ON THE LIMIT DISTRIBUTIONS OF CONTINUOUS-STATE BRANCHING 

PROCESSES WITH IMMIGRATION 



MARTIN KELLER-RESSEL AND ALEKSANDAR MIJATOVIC 

Abstract. We consider the class of continuous-state b ranchins processes with immigration (CBI- 
processes), introduced bv lKawazu and Watanabd [l97l| and their limit distributions as time tends 



to infinity. We determine the Levy-Khintchine triplet of the limit distribution and give an explicit 
description in terms of the characteristic triplet of the Levy subordinator and the scale function of 
the spectrally positive Levy process, which describe the immigration resp. branching mechanism 
of the CBI-process. This representation allows us to describe the support of the limit distribution 
and characterise its absolute continuity and asymptotic behavior at the boundary of the support, 
generalizing several known results on self-decomposable distributions. 



1. Introduction 

Continuous-state branching p rocesses with immigration (CBI-processes) have been introduced 



Kawazu and Watanabe 



197111 as sca l ing li mits of discrete single-type branching processes with 



Kawazu and Watanabe 



197l[ | the authors show that in general a CBI-process has a 



by 

immigration. In 

representation in terms of the Laplace exponents F{u) and R{u) of two independent Levy processes: 
a Levy subordinator and a spectrally positive Levy process X^, which can b e interpreted as 
immig ration and branching mechanism of the CBI-process respectively. Following 



Dawson and Li 



2006l | a CBI-process can in fact be represented as the unique strong solution of a non-linear SDE 



driven by X^ and X^. 

For discrete branching; processes with immigration, limit distributions have been studied already 



m 



Heathcote 



1965 



196611 and so me results on the existence of limit distributions of a CBI- process 







Pinskyl 19721 ] . albeit witho ut proofs. Re cently, proofs for the results of iPinsky 



201 1[ | . The main result of lPinskyI |l972l | states that under an integral condition 



were published bj^ 
have appeared in 

on the ratio F[u)/R{u) a limit distribution exists and can be described in terms of its Laplace 
expon ent (cf. Theorem I2.6p . The contribution of this article is to build on the results of Ipinskv 



1972l | in order to give a finer description of the limit distribution: We show that it is infinitely 
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divisible, give a representation of its Levy-Khintchine triplet (Theorem 13. ip and then use this 
new representation to obtain results on smoothness, support and other properties of the limit 
distribution. From this main result, several other representations of the Levy-Khintchine triplet 
are then derived. The most concise representation is given by equation (j3.17p . which states that the 
Levy measure of the limit distribution has a density of the form x i— > k{x) /x and the corresponding 
fc-function k : (0, oo) — )• M^o is given by the formula 



k = -Aj^pW, 

where A-^f is the generator of the modified Levy subordinator and W is the scale function 
that corresponds to the spectrally one-sided Levy process see Remark 13.41 for the precise 

statement of this factorization). In Section [J] we derive further properties of the limit distribution. 
In particular, we characterize the support of the limit distribution, show its absolute continuity 
and describe its boundary behavior at the left endpoint of the support. Furthermore we prove that 
the class of limit distributions of CBI-processes is strictly larger than the class of self-decomposable 
distributions on M^o and is strictly contained in the class of all infinitely divisible distributions on 

Most of our results can also be regarded a s extensions of known res u lts on limit distributions o f 
Qrnstein-U hlenbeck-type (OU-type), see e.g. Jurek and Vervaat 1983 |. [sato and Yamazato 1984 1. 



Sato 



1999(1, to the class of CBI-processes. The knowledge of the Levy-Khintchine triplet for station- 



ary distrib utions of OU-type processes has been applied to statistical estimation of the underlying 



Masuda and Yoshida 



process m 

used for extensions of this methodology to CBI-processes. 



20051]). We suggest that in further research our results may be 



2. Preliminaries 



2.1. Continuous-State Branching Processes with Immigration. Let X b e a continuous 



state branching process with immigration. Following iKawazu and Watanabd [I97l[, such a process 
is defined as a stochastically continuous Markov process with state space [0, oo], whose Laplace 
exponent is affine in the state variable, i.e. there exist functions (p{t, u) and ili{t, u) such that 



(2.1) 



logE^ [e 



(/>(t, u) + xV'(t, u) 



for all i > 0, li > 0, a; G 



where, as usual for the theory of Markov processes, denotes expectation, conditional on Xq = x. 
Since we are interested in the limit behavior of the process X as t ']" oo, we further assume that X is 
conservative, i.e. that Xt is a prope r random variable for each t > with state space M3.0 := [0, 00). 
The following theorem is proved in lKawazu and Watanabd 19711 ]. 
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Theorem 2.1 (jKawazu and Watanabd 197ll |). Let {Xt)t>o be a conservative CBI-process. Then 
the functions (j){t, u) and ip{t, u) in (|2.ip are differentiahle in t with derivatives 



(2.2) 



F{u)= 



t=o 



R{u) = ^^^{t,u) 



t=0 



and F , R are of Levy-Khintchine form 

(2.3) F{u) = bu- [ (e-"« - l) m{d(), 

(2.4) R{u) = -au^ + /3n - / (e-< - 1 + </(o,i] iO) K^O, 

J(0,oo) ^ ^ 

where a, 6 € ^^o, /3 G M, /(o,i] is the indicator function of the interval (0, 1] and m, fi are Levy 
measures on (0,oo), with m satisfying f(^Q^-^ {x A l)m(dx) < oo, and R satisfying 



(2.5) 



1 



-ds = oo where R*{u) = max(i?(ti), 0). 



'0+ R*{s) 

Moreover (j){t,u), ip{t,u) take values in M^q and satisfy the ordinary differential equations 

d 



(t){t,u) = F {^l:{t,u)) 



0(0, ^z) = 0, 



dt 

^^ij{t,u) = R{^{t,u)) , ^(0,n)=«. 



(2.6) 



Remark 2.2. The equations (j2.6p are often called generalized Riccati equations, since they are 
classical Riccati differential equations, when m = fi = 0. 

We call (F, R) the function al characteristics of the CBI-process X. Furthermore the article of 



Kawazu and Watanabe 



19711 ] contains the following converse result: for any functions F and R 



defined by (j2.3p and (j2.4p respectively, which satisfy condition (j2.5p and the restrictions on the 
parameters a,b, f3 and the Levy measure m stated in Theorem 12.11 there exists a unique con- 
servative CBI-process with functional characteristics {F, R) . In this sense the pair (F, R) truly 
characterizes the process X. Clearly, F{u) is the Laplace exponent of a Levy subordinator X^ , 
and R{u) is the Laplace exponent of a Levy process X^ without negative jumps. Thus, we also 
have a one-to-one correspondence between (conservative) CBI-processes and pairs of Levy pro- 
cesses {X^,X^), of which the first is a subordinator, and the second a process without negative 
jumps that satisfies condition ()2.5p . In the case of a CBI-process without immigration (i.e. a 
CB-process), which correspon ds to F = 0, a pathwise transformation of X^ to X and vice versa 
was given bv iLampertil 1967l |. and is often referred to as 'Lamperti transform'. Recently, a path- 
wise correspondence between the pair (X^ , X^) and the CBI-process X has been constructed by 
Caballero. Perez Garmendia. and Uribe Bravd 201C]. 



-'^The notation J denotes an integral over an arbitrarily small right neighborhood of 0. 
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The following properties of F{u) and R{u) can be easily derived from the representations ()2.3p 
and (j2.4p and the parameter conditions stated in Theorem 12. li 

Lemma 2.3. The functions F{u) and R{u) are concave and continuous on Mj>o and infinitely 
differentiable in (0,oo). At u = they satisfy F{0) = R{0) = 0, the right derivatives -F^(O) and 
-R^(O) exist in (— cx),+oo] and satisfy F!^{0) = liuiu^oF'^u) and R'^{0) = lim„j^o ^'(''^)- 



We will also need the following result, which can be found e.g. in iKvprianoul 20061 . Ch. 8.1] 

Lemma 2.4. For the function R{u) exactly one of the following holds: 

(i) R'^{0) > and there exists a uq > such that R{uq) = 0; 

(ii) R = 0; 

(Hi) i?'+(0) < and R{u) < for all u>0. 

Remark 2.5. In case ^ R{u) is called a supercritical branching mechanism, while case can be 
further distinguished into critical (i?^(0) = 0) and subcritical branching {R'_^{0) < 0). 

In what follows we will be interested in the limit distribution and the invariant distribution 
of {Xt)t>o- We write Ptf{x) = E^[f{Xt)] for all x G R^o and denote by {Pt)t>o the transition 
semigroup associated to the Markov process X. We say that L is the limit distribution of the 
process X = (Xt)t>o if Xt converges in distribution to L under all P^' for any starting value 
x € M^o of X. We call L an invariant (or stationary) distribution of X = {Xt)t>o, if 



[0,oo) 

for any t >0 and bounded measurable / : 
Lby 

l{u) = - log 



Ptf{x)dL{x) 



/(x)dL(x) 



Finally we denote the Laplace exponent of 



[0,oo) 



'dL{x) 



{u > 0). 



2.2. Limit Distributions of CBI-Processes. Theorem 12.61 and Corollary 12.81 concern the ex- 
istence of a limi t distribution of a CBI-process and have been anno unced in a similar form but 
without proof in Pinskv 1972]. A proof has recently appeared inQ [20 111 . Thm. 3.20, Cor 3.21]; 
the only difference to the result given here is that we drop a mild moment condition assumed in 



12011 



, Eq.(3.1)f] and that we include stationary distributions in the statement of our result. 
S ome weaker results on the existence of a limit distribution of a CBI-process have also appeared 
in 



Keller-Ressel and Steiner 



20081 ]. We give a self-contained proof of the theorem and its corollary 



in the appendix of the article. 

Theorem 2.6 (|pinskvl [iOtJ^Q |201i| ). Let {Xt)t>o be a CBFprocess on . 
statements are equivalent: 

(a) {Xt)t>Q converges to a limit distribution L as t 



Then the following 



00: 
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(b) {Xt)t>o has the unique invariant distribution L; 

(c) It holds that R'+iO) < and 

(2.7) - / -^ds < oo 



R{s) 
for some u > 0. 
Moreover the limit distribution L has the following properties: 
(i) L is infinitely divisible; 

(a) the Laplace exponent l{u) = — log /[q oo) e~"^dL(2;) of L is given by 

(2.8) l{u) = - / -y ds (n > 0) . 

Jo R[s} 



Remark 2.7. Note that the existence of the right derivative i?^ at that appears in statement (jcj) 
is guaranteed by Lemma 12.31 

Corollary 2.8. // R'j^{0) < then the integral condition (j2.7p is equivalent to the log-moment 
condition 

(2.9) / logCm(dO<oo. 



2.3. Results on Ornstein-Uhlenbeck-type Processes. A subclass of CBI-processes, whose 
limit distributions have been studied extensively in the literature is the class of M^o-valued Ornstein- 
Uhlenbeck-type (OU-type) processes. We briefly discuss some of the known results on OU-type 
processes, that will be generalized by our results in the next section. Let A > and Z he a Levy 
subordinator with drift b € M^q ai^d Levy measure 7?T,(d.^). An M^q" valued OU-type process X is 
the strong solution of the SDE 



(2.10) dXt = -XXtdt + dZt, Xq € 



which is given by Xt = Xqc"'*'* -|- e^^'^~^^dZs. This is the classical Ornstein-Uhlenbeck process, 
where the Brownian motion has been replaced by an increasing Levy process. It follows from 
elementary calculations that an M^Q-valued OU-type process is a CBI-process with R{u) = —Xu. 
In terms of the two Levy processes X^ , X^, this corresponds to the case that X^ = Z, and X^ 



is the degenerate Levy process Xj^ = —Xt. For OU-type proce sses analogues of 



leorem 12.61 and 



197ll ]. 



the log-moment condition of Corollary 12.81 already appeared in iQinlar and Pinskv 

An interesting characterization of the limit distributions of OU-type processes is given in terms 
of self-decomposability: Recall that a random variable Y has a self- decomposable distribution if for 
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every c G [0, 1] there exists a random variable Yc, independent of Y, such that 



(2.11) 



Y = cY + Yc. 



Self-decomposable distributions are a subclass of infinitely divisible distributions, and exhibit in 
many aspects an increased degree of regularity. It is known f or exa r nple, that every non-degenerate 
self-decomposab le d istribution is absolutely continuous (cf. ISatd |1999l . 27.8]) and unimodal (cf. 



Yamazato 



19781 1 or 



Sato 



19991 . Chapter 53]), neither of which holds for general infinitely divisi- 



ble distributions. As we are working with non-negative processes, we focus on self-decomposable 
distributions on the half- line M^Qj and we denote this class by SD+. The connection to OU-type 
processes is made by the following result. 



Theorem 2.9 (jJurek and VervaatI [19831 ]. ISato and Yamazatd [19841 ]). Let X be an OU-type pro- 
cess on R^o onii suppose that m(d^) satisfies the log-moment condition J^^-j^ log^m(d^) < oo. 
Then X converges to a limit distribution L which is self-decomposable. Conversely, for every 
self-decomposable distribution L with support M^q there exists a unique subordinator Z with drift 
b G M^o '^iT'd a Levy measure m(d^), satisfying J^^-^ log^?7i(d^) < oo, such that L is obtained as 
the limit distribution of the corresponding OU-type process. 



Since a self-decomposable distribution is infinitely divisible, its Laplace exponent has a Levy- 
K hintchine d ecomposition. The following characterization is due to Paul Levy and can be found 
m Isatol |l999l . Cor. 15.11]: an infinitely divisible distribution L on M^o is self-decomposable, if and 
only if its Laplace exponent is of the form 



(2.12) 



log / e 

'[0,oo) 



'^dL{x) = 'ju 



where 7 > and A; is a decreasing function on 
Levy subordinator Z by 



1) M^dx, 



The parameters 7 and k are related to the 



(2.13) 



and 



k{x) 



—m{x, 00) 
A 



for X > 0. 



Following ISatd [19991 ] we call k the fc-function of the self-decomposable distribution L. Many prop- 
erties of L, such as smoothness of its density, can be characterized through k. In fact, several sub- 
classes of SD_|_ have been defined, based on more restrictive assumptions on k. For example, the class 
of self-decomposable distributions whose fc-function is completely monot one, is known as the Thorin 
class, and arises in the study of mixtures of Gamma distributions; see iJames. Rovnette. and Yor 



20081 ] for an excellent survey. In our main result. Theorem 13.11 we give analogues of the formulas 



(j2.12p and (j2.13p for the limit distribution of a CBI-process. As it turns out, a representation as 
in (|2.12p still holds, with the class of decreasing fc- functions replaced by a more general family. 
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However, we do not obtain a structural characterization of the CBI hmit distributions that re- 
places self-decomposability. Identifying such a structural condition (if there is any) constitutes an 
interesting question that is left open by our results. 

3. Levy-Khintchine Decomposition of the Limit Distribution 

Let and be the Levy processes that correspond to the Laplace exponents F and R 
given in (j2.3p and (j2.4p respectively. As remarked in Section 12.11 X^ is a subordinator and X^ 
is a Levy process with no negative jumps. From Theorem 12.61 it follows that whenever a limit 
distribution exists, then E [Xf] = R+{0) < and i? ^ 0, such that X^ is not a subordinator. 



but a true spectrally positive Levy process in the sense of iBertoinI 1996]. The fluctuation theory of 
spectrally one-sided Levy processes has been studied extensively. The convention used in much of 
the literature is to study a spectrally negative process. In our setting such a process is given by the 

R( u) for ^ ^ 0. A central result in 



dual X^ = —X^ and its Laplace exponent is logE e^^' 
the fluctuation theory of spectrally one-sided Levy processes (see [Bertoinl . Il996l . Thm 8, Ch VII]) 
states that for each function R of the form ()2.4p there exists a unique function : M — ?• [0,oo), 
known as the scale function of X^, which is increasing and continuous on the interval [0, oo) with 
Laplace transform 

1 



(3.1) 



f 

Jo 



e-"^W^(x)dx 



'i?(n)' 



for n > 0. 



and identically zero on the negative half-line {W{x) 
implies that 



for all X < 0). Note that this equality 



(3.2) 



W{x) = o{e''^) as x — )• oo for any e > 0, 



i.e. that W has sub-exponential growth, a fact that will be needed subsequently. Furthermore the 
scale function W has the representation 
W{x) ( rv 



(3.3) 



W{y) 



exp 



n{£'> z) dz 



for any 



< X <y, 



where n is the Ito excursion measure on the set 



(3.4) £ = {ee D(R) : 3Ce G (0, oo] s.t. e{t) = if Ce < t < oo, e(0) > 0, £{t) > OVt e (0, (s)} , 

with L'(R) the Skorokhod space. The measure n is the intensity measure of the Poisson point 
process of excursions from the supremui n of X^ and j e > z} C £ denotes the set of excursions of 



height e = sup^<^^ e(t) at least z > (see lBertoinI [19961 ] for details on the Ito excursion theory in the 
context of Levy processes). The representation ()3.3p implies that, on the interval (0, oo), the scale 
function W is strictly positive, absolutely continuous, log-concave with right- and left-derivative 
given by W!^{x) = n{e > x) W{x) and WLix) = n{e > x) W{x) respectively. Furthermore at x = 
the right-derivative W+(0) exists in [0, cxo]. 
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Using the scale function W associated to we can formulate our main result on the Levy- 
Khintchine decomposition of the limit distribution of a CBI-process. 



Theorem 3.1. Let X be a CBI-process with functional characteristics {F,R) given in (j2.3p 
and (|2.4|) and assume that X converges to a limit distribution L. Let W be the scale function 
associated to the dual X^ of the spectrally positive Levy process X^ and let (6, m) be the drift and 
Levy measure of the subordinator X^ . Then L is infinitely divisible, and its Laplace exponent has 
the Levy-Khintchine decomposition 



oo 



(3.5) - log / e-"ML(x) = uj- [ (e"^" - l) ^ dx, 

Jo J(0,oo) X 

where 7 > and k : (0, 00) are given by 

(3.6) 7 = bW{0), 

(3.7) k{x) = bW'^{x)+ f [W{x)-W{x-i)]m{di). 

J(0,oo) 

Remark 3.2. If X is an OU-type process, then R is of the form R{u) = —An with A > 0. Since 
the scale function in this case takes the form W{x) = ;x-^[o,oo)(^); equations (|3.6p and (13. 7p reduce 
to 7 = and k[x) = ■^m(x,oo). This is precisely the known result for the M^o-valued OU-type 
process stated in (|2.13p . 

Before proving this result, we state a corollary that connects the limit distribution in Theorem 13. II 
with the excursion measure n associated to the Poisson point process of excursions away from the 
supremum of the dual of the branching mechanism X^. To state it, we introduce the effective drift 
Ao of X^, which is defined as 



(3.8) 



/{o 1] ^^('^?) ~ Z^' bounded variation. 



+00, if X^ has unbounded variation. 

Note that Aq > must hold if R'+{Q) < 0. 

Corollary 3.3. Let the assumptions of Theorem \3. 1\ hold. Then 

(3.9) 7 ~ 1~ 

Ao 

(3.10) k{x) =W{x) ^bn{e> x) + J (^l - l{^<r,}exp J n{e > z)dz^^ m{dC)^ 

where Aq is the effective drift of X^, n is the ltd excursion measure corresponding to the Poisson 
point process of excursions from the supremum of X^ . 
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Proof of Theorem \3.1[ For every t > 0, the distribution of Xt is infinitely divisible and supported 
on M^o- Hence the same is true of the limit distribution L. The Laplace exponent of L can by 
Theorem 12.61 be expressed as 

(3.11) -log/ e-"^dL(x) = -/ ^^ds = du- [ (e""^ - l) i/(dx), u>0, 

J[0,oo) Jo J{0,oo) 

for d > and some Levy measure i^(dx) satisfying J^^ xi'{dx) < oo. Moreover, it is clear from 
Theorem [M] that i?'+(0) < and ^ 0. Thus, Lemma [13] implies that the quotient F/R is 
continuous at any n > 0. Since the elementary inequality \e~^^ — l\/h < x holds for all x,h > 0, 
the dominated convergence theorem and the fundamental theorem of calculus applied to (j3.1ip 
yield the identity 

(3.12) -^^7^ = d+ [ e-"^xi/(da;) for ah u > 0. 



R{u) 7(0,00) 

Any twice-differentiable function / that tends to zero as — t- oo, i.e. / G Cq{M.), is in the 
domain of the generator Axf of the subordinator and the following formula holds 

(3.13) Axf fix) = bf\x)+ [ [/(x + e)-/(x)] m(dC) for xeR. 

J(0,oo) 

Fix u > and let /„ G C'o (^) ^ function that satisfies fu{x) = e~"^' for all x > 0. Applying ()3.13p 
to fu yields Ax^fuix) = F{u)fu{x) for all ti > 0, x > 0. Multiplying by —W{x) and integrating 
from to oo gives the following identity for all n > 0: 

/•oo 

(3.14) - / W{x)AxFfu{x)dx 

Jo 

= b ue-"^VI/(x)dx- / W{x) / e"<^+0_e-ux ^^(d^)^^;. 

Jo Jo J{0,oq) L -I 

Note that AxFfuix) ~ e~"^ for x ^ oo, which guarantees that the integrals are finite in light of 
eq. (|3.2p . Since W is increasing and absolutely continuous, integration by parts gives 

f'OQ rOQ 

(3.15) - / W{x)AxFfu{x) dx = bW{0) + b e-"^H^+(x) dx 

Jo Jo 

+ / e-"^ / [W{x)-W{x-0]m{dC)dx, 

Jo J{0,oq) 

for all n > 0. The second integral on the right-hand side of p.lSp is a consequence of the following 
steps: (i) note that the corresponding integrand in (|3.14p does not change sign on the domain of 
integration, (ii) approximate the Levy measure m by a sequence of measures (m„)„gN with finite 
mass, (iii) apply Fubini's theorem to obtain the formula for each m„, (iv) take the limit by applying 
the monotone convergence theorem. 
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On the other hand, combining the identity Axpfuix) = F{u)e for all n > 0, x G with 
311 and (fSHT) yields 

Fiu) 



(3.16) 



W{x)AxFfu{x)dx 



VF(x)F(n)e~"^dx 



d 



^xv{(lx), 



10 JO R{u) 7(0,oo) 

which in turn must equal the right hand side of (|3.15p . In the limit as n ^ 00, the equality of the 
expressions in ()3.15p and ()3.16p yields d = bW{0). Subtracting this term we arrive at the equality 



/ e"'^^ bW[{x)+ / [W{x)-W{x-^)]m{d^) dx = 

Jo J{0,oo) J(i 



'(0,00) 



e '^^xi^{dx). 



Both sides are Laplace transforms of Borel measures on (0, 00), and we conclude from the equality 
of transforms the equality of the measures 



dx = xi'{dx) 



bW+{x) + / [W{x) - W{x - C)] m{d^) 
(0,00) 



for all x > 0. In particular it follows that i/(dx) has a density k{x)/x with respect to the Lebesgue 
measure, and that k{x) is given by 



□ 



k{x) = bW'^{x) + W{x)m{x, 00) + / [W{x) - W{x - i)] 

J{0,x] 



almost everywhere, which concludes the proof. 



Proof of Corollary Following iKvprianoul 20061 . Lemma 8.6], VF4.(0) > if and only if has 
finite variation, and is equal to I/Aq in this case, with Aq defined by (j3.8p . If X^ has infinite 
variation, then 1^(0) = and Aq = 00, and hence equation (j3.9p holds. Substituting the represen- 
tation (j3.3p of W in terms of the measure n into (j3.7p yields the second equation (j3.10p . □ 



Remark 3.4. The formula for the /c-function in (j3.7p looks very much like the Feller generator of 
the subordinator X^ applied to the scale function W of X^. However, the Feller generator is only 
defined on a subset (i.e. its domain) of the Banach space of continuous functions that tend to at 
infinity, Co(K). Any function in the domain of the generator of X^ must be in Cn(M) and differen - 



tiable ; sufficient conditions for the differentiability of W are given in lChan. Kvprianou. and Savov 
201l[. However, the scale function W, which is non-decreasing on M^O) is not in Co(M) and thus 



never in the domain of the Feller generator of X^ . To remedy this problem, consider that by (j3.2 
and p.3p . both W and W!^ are elements of the weighted Li-space defined by 



L?(0,«)) :=<j/GL'r(0,oo) 



\f{x)\h{x)dx < 00 



where h : (0, 00) (0, 00) is a continuous bounded function with lim2,.|o h{x) = and h{x) 



~ e 



as X 
dl 



00 for some c > 0. The semigroup {Pt)t>o of the Markov process X- 



F,x 



{Xf<x} 



+ 



{X[>x} 



(i.e. the dual, started at x > 0, of X , sent to a killing state d upon the first passage 
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E 



f 



F,x 



E 



/ {X - X[) l^^.^^y 



for each 



into (—00,0)) acts on LJ'(0, 00) by Ptf{x) 
f € L^(0,oo) where we take f{d) = 0. It can be shown that the -L^(0, 00) -semigroup (-fj;)t>o is 
strongly continuous with a generator A-^p and, if -R'|_(0) < 0, then the scale function W associated 
to R is in the domain of ^j^^. Furthermore the /c-function k in Theorem 13.11 can be written as 

(3.17) k = -Aj^fW. 

The proof of these facts is straightforward but technical and rather lengthy and hence omitted. 



4. Further Properties of the Limit Distribution 

As discussed in Section [Z3l self-decomposable distributions, which arise as limit distributions of 
M^o-valued OU-type processes are in many aspects more regular than general infinitely divisible 
distributions. For self-decomposable distributions precise results are known about their support, 
absolute continuity and behavior at the boundary of their support. Using the notation of Sec- 
tion 12.31 and excluding the degenerate case of a distribution concentrated in a single point, the 
following holds true when L is self-decomposable: 

(i) the support of L is [b/X,oo); 

(ii) the distribution of L is absolutely continuous; 

(iii) the asymptotic behavior of the density of L at 6/A is determined by c = limj-^o k{x). 
We refer the reader to 



Sato 



19991 . Theorems 15.10, 24.10, 27.13 and 53.6]. The goal in this section 
is to show analogous results for the limit distributions L arising from general CBI-processes, i.e. to 
characterize the support, the continuity properties and the asymptotic behavior at the boundary 
of the support of L. We start by isolating the degenerate cases. A Levy process, such as or 
X^, is called degenerate, if is deterministic, or equivalently if its Laplace exponent is of the form 
u Xu for some A € M. For a CBLprocess X we draw a finer distinction. 

Definition 4.1. A CBI-process X is degenerate of the first kind, if it is deterministic for all starting 
values Xq = 2; € M^jq. X is degenerate of the second kind, if it is deterministic when started at 
Xo = 0. 

Remark 4.2. Clearly degeneracy of the first kind implies degeneracy of the second kind. From 
Theorem 12.11 the following can be easily deduced: a CBLprocess is degenerate of the first kind if 
and only if both X^ and X^ are degenerate. In this case F{u) = bu and R{u) = f3u, and X is the 
deterministic process given by 

(4.1) X, = Xoe^' + ^ (e^* - 1 
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A CBI-process X is degenerate of the second kind, but not of the first, if and only if = and 
X^ is non-degenerate. In this case it is a CB-process, i.e. a continuous-state branching process 
without immigration. 

The following proposition describes the support of the limit L in the degenerate cases. 

Proposition 4.3. Let X be a CBI-process, let L he its limit distribution and let k be the function 
defined in Theorem \3.1[ If X is degenerate of the first kind, then suppL = {—b/j3}. If X is 
degenerate of the second kind but not of the first kind, then suppL = {0}. Moreover, the following 
statements are equivalent: 

(a) the support of L is concentrated at a single point; 

(b) X is degenerate (of either kind); 

(c) there exists a sequence Xj J, such that k{xi) = for all i € N; 

(d) k{x) = for all x > 0. 

Proof. Suppose that X is degenerate of the first kind. Then the limit L is concentrated at —b//3 
by (|4.ip . Suppose next, that X is degenerate of the second, but not the first kind. Then F = 
and by Theorem 13.11 the Laplace exponent of L is 0. It follows that L is concentrated at in this 
case. 

We proceed to show the second part of the proposition. It is obvious that ([d]) implies (jcj). To 
show that (jcj) implies (|b|), note that the inequality 

k{x) > W{x) {hnie > x) -|- m(x, oo)) 

holds for all X > by equation (I3.10p . Since W{x) > for any x > 0, assumption (jcj) implies that 

bn{e > Xi) + m{xi, oo) = for all Xj, « G N. 

We can conclude that m = and hence Xf = bt for all t > 0. Furthermore we see that either 
6 = or n = 0. If 6 = 0, then F = and hence, by Remark 14. 2j X is a degenerate CBI-process of 
the second kind. On the other hand, if the Ito excursion measure n is zero, then the representation 
in (jS.Sp implies that the scale function W is constant. In this case it follows from (j3.ip that 
R{u) = (3u for some /3 < 0, or equivalently that X^ = fit for all t > and hence that X is 
degenerate of the first kind. 

The fact that (Jb]) implies (jaj) follows from the first part of the proposition. It remains to show 
that (jaj) implies (jdj); this is a consequence of the fact that L is infinitely divisible with support in 
M^O) and that the support of an infinitely divisibl e distributio n in M^o is concentrated at a single 



point if and only if its Levy measure is trivial (cf. ISatd [19991 . Thm. 24.3, Cor. 24.4]). □ 



The next result describes the support of the limit L in the non-degenerate case. 
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Proposition 4.4. Let X be a non- degenerate CBI-process and let L he its limit distribution. Then 

suppL = [6/Ao, oo), 

where Aq is the effective drift of X^, defined in (|3.8p . In particular suppL = Mj>o if and only if 
b = or the paths of X^ have infinite variation. 

Proof. From Proposition 14.30 we know that there is some 6 > 0, such that the /c-function of L is 
non-zero on (0, 6). For any h G (0, 5), define Lh as the infinitely divisible distribution with Laplace 
exponent J^^ (e~^" — 1) ^^dx. Each Lf^ is a compound Poisson distribution, with Levy measure 

I'hidC) = ^^l(h,oo)(2;)- Since k is non-zero on (0,5) 

(4.2) {h,6) C suppi/h C (/i, oo). 

From Theorem 13.11 we deduce that as h the distributions converge to L{j + .), i.e. to L 
shifted to the left by 7. For the supports, this implies that 



(4.3) suppL = {7} -|- limsuppL/i 

where the limit denotes an increasing union of sets and '-|-' denotes pointwise addition of sets. 
Using ()4.2p and the fact that is a compound Poisson distribution it follows that 



{0}U [J{nh,n6) C suppL/j C {0} U (/i, 00) 



n=l 



bv lSatol [19991 . Thm. 24.5]. Let h 10 and apply ([43]) to obtain 

00 

[Jb,7 + n6] C suppL C [7,00), 

n=l 

and we conclude that suppL = [7, 00). By Corollarv 13.31 = 6/Ao, which completes the proof. □ 



Proposition 4.5. Let X be a CBI-process and let Aq be as in (|3.8p . Then the limit distribution L 
is either absolutely continuous on M^o or absolutely continuous on M^o \ {V-^o} with an atom at 
{6/Ao}, according to whether 

(4.4) / ^-^dx = 00 or / —-^dx < 00. 

Jo X Jo X 

Proof. If X is degenerate, then the assertion follows immediately from Proposition 14. 3i In this case 
k{x) = for all x > 0, the integral in (j4.4p is always finite and the distribution of L consists of a 
single atom at 6/Ao. 

It remains to treat the non-degenerate case. Assume first that the integral in (14. 4p takes a finite 
value. Then also the total mass ^{0, 00) of the Levy measure i/(dx) = ^^dx is finite, and L — 'j 
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has compound Poisson distribution. By ISatd 19991 . Rem. 27.3] this imphes that for any Borel-set 

A C Ms^o 



oo 



(4.5) / dL(a;) = e-*'^(0'~)^^z.*^(^), 

where v*^{(ix) is the j-th convolution power of v, and it is understood that v*^ is the Dirac measure 
at 0. Since z/(dx) is absolutely continuous - it has density - also the convolution powers i'*^{dx) 
are absolutely continuous for j > 1. The first summand v*^ however has an atom at 0. It follows 
by ()4.5p that L — 7 is absolutely continuous on (0, 00) with an atom at 0, and we have shown the 
claim for the case ^^dx < 00. 

Assume that ^^ dx = 00. Then the Levy measure v{dx) = ^^dx of L has infinite total 



mass, and ISatd 1999I . Thm. 27.7] implies that L has a distribution that is absolutely continuous, 
which completes the proof. □ 

So far, we know that the left endpoint of the support of L is 7 = 6/Ao, and that the distribution 
of L may or may not have an atom at this point. In case that there is no atom, the following 
proposition yields an even finer description of the behavior of the distribution close to 7. 

Proposition 4.6. Let X he a CBI-process satisfying the assumptions of Theorem \3.1[ and let L 

he its limit distribution. Suppose that c = lim^jio fc(x) is in (0, 00), and define 

(4.6) Kix) = exp(^l\c-kiy))^y 
Then K{x) is slowly varying at and L satisfies 

(4.7) L(x) ~ -^(x-7)^-ii^(x-7) as x 1 7, 

r(c) 

where 7 = 6/Ao and 

/ /■! r°° dx f°° dx\ 

Proof. Note that the inequality c > and Proposition 14. 5 1 imply that L is abso l utely continuous. Its 



support is by Proposition 14.41 equal to [7, 00) and ^(7) = 0. The proof of lSatd jl999l . Theorem 53.6] 



for self-decomposable distributions can now be applied without change. □ 

Recall that ID_|_ and SD_|_ denote the classes of infinitely divisible and self-decomposable distri- 
butions on M^o respectively. Let CLIM be the class of distributions on Mj,o that arise as limit 
distributions of CBI-processes. 

Proposition 4.7. The class CLIM is contained strictly between the self-decomposable and the 
infinitely divisible distributions on M^O; ^-e. 

SD . C CLIM C ID . . 
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Proof. The inclusion CLIM C ID+ follows from Theorem 12. 6| and the inclusion SD+ C CLIM from 
Theorem 12.91 and the fact that each M^o-valued OU-type process (see (|2.10p ) is a CBI-process with 
R'-^{0) = —X < 0. The strictness of the inclusions can be deduced from the following facts: 

• all distributions in SD+ are either degenerate or absolutely continuous (cf. ISatd 1999I . 
Thorem 27.13]); 

• all distributions in CLIM are absolutely continuous on M^jq \ {b/Xo}, but some concentrate 
non-zero mass at {6/Ao} (cf. Propositions 14.41 and [ 

• the class ID+ contains singular distributions (cf. 



Sato 



19991 . Theorem 27.19]). 



□ 

For a more direct proof of the fact that SD+ is strictly included in CLIM we exhibit an example 
of a distribution that is in CLIM but not in SD+: 

Example 4.8 (CBI-process with non self-decomposable limit distribution). In this example we con- 
sider the class of CBI-processes X given by a general subordinator and spectrally positive 
process X^ equal to a Brownian motion with strictly negative drift. The Laplace exponent of X^ 
is R{u) = —ail} -|- /3ti with a > 0, /? < 0. It is easy to check using (|3.1|) that the scale function of 
the dual X^ and its derivative are 

(4.8) W{x) = [exp (x/3/a) - 1] //3 and W'{x) = exp {xf3/a) /a. 

Theorem 13.11 implies that the characteristics of the limit distribution L are given by 7 = and 



(4.9) k{x) = e^^/" 




m(x, 00) -|- 



1 



m(d.^) 



m(x, oo)//3, 



where b G M^o is the drift and m the Levy measure of the subordinator X^ . Assuming in addition 
that X^ is a compound Poisson process with exponential jumps and setting parameters equal to 



m{x, 00) 

formula (|4.9|) reduces to k{x) = 2(e~^ 



0, a = 1/2, (3 



-h 



-2x\ 



Since this fc-function is not decreasing, the corre- 



sponding distribution L, which is in CLIM, cannot be in SD-|_. 

Proposition (j4.9|) gives sufficient conditions for a distribution in CLIM to be self-decomposable. 



Proposition 4.9. Let X be a CBI-process and let L be its limit distribution. Each of the following 
conditions is sufficient for L to be self-decomposable: 

(a) /.f = and a = 0, 

(b) 11 = and m = 0, 

(c) m = and W is concave on (0,oo). 

Conversely, if m = and L is self- decomposable, then W must be concave on (0,oo). 
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Remark 4.10. The monotonicity of the derivative of the scale function, which ari ses in Proposi - 
tion 14.91 also plays a role in other applications of scale functi ons (e.g. control theory iLoeffenl 2008l |: 



conjugate Bernstein functions and one-sided Levy processes iKvprianou and Rivero 



20081 1 1. 



Proof. The first two conditions are rather trivial. In the first case X is an QU-type process, and self- 
decom posability follows from the classical results of 



Jurek and Vervaat 



19831 1 



Sato and Yamazato 



1984i | that we state as Theorem 12.91 In the second case X has no jumps, and hence is a Feller 



diffusion. This process is well-studied, and its limit distribution is known explicitly. It is a shifted 
gamma distribution, which is always self-decomposable. It remains to show (c) and the converse 
assertion. Assume that m = 0. By Theorem 13. II we have k{x) = bW!^{x) in this case. An infinitely 
divisible distribution is self-decomposable if and only if it can be written as in (j2.12p with decreasing 
/c-function. Clearly k is decreasing if and only if WL is, or equivalently if W is concave on M^q. □ 



Appendix A. Additional Proofs for Section [2] 

Proof of Theorem \2M We first show that (|cj) is equivalent to (jaj) and that L has to satisfy (ji]) 
and dni). Consider the three alternatives for the behavior of R{u) that are outlined in Lemma l2.4i 
Through the Riccati equations (|2.6p they imply the following behavior of tp{t,u): If -R+(0) > 
then limt^oo ^(i, u) = uq for all t,u > 0, if R = then ip{t, u) = u for all u > 0, and if ii+(0) < 
but i? ^ then lim^^oo "^{i-, u) = 0. Moreover, 

POD 

(A.l) lim - logE [e"''^*] = lim {(j){t, u) + xTp{t, u)) = / F{7p{r, u)) dr + x ■ lim tpit, u). 

t\oo t\oo Jq t^oo 

We see that if -R+(0) > or = the right-hand side diverges for n > 0, and hence no limit 
distribution exists in these cases. In case that -R+(0) < and R^ 0, the transformation s = ijj{r, u) 
yields that 

(A.2) lim-logE fe-"^*] = / ^^ds. 

This integral is finite, if and only if condition (j2.7|) holds. If it is finite then Levy's continuity theorem 
for Laplace transforms guarantees the existence of, and convergence to, the limit distribution L 
with Laplace exponent given by (j2.8p . It is also clear that L must be infinitely divisible, since it is 
the limit of infinitely divisible distributions. If on the other hand the integral in ()A.2p is infinite 
for some u G M^Qj then there is no pointwise convergence of Laplace transforms, and hence also no 
weak convergence oi Xt as t ^ oo. 

To complete the proof it remains to show that any limit distribution is also invariant and vice 
versa, i.e. that (jb]) is equivalent to (jaj). Assume that L is an invariant distribution of {Xt)t>o, 
and has Laplace exponent l{u) = — log e~"^d-L(x). Denote fu{x) = e~"^ and note that the 
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invariance of L implies 

(A.3) / fu{x)dL{x) = I PtUx)dL{x) = e--^^*'") / e-^^(*'")dZ(x) 

J [0,00) J [0,00) J [0,00) 

for all t, n > 0. This can be rewritten as l{u) = (j){t, u) + l{ip{t, u)). Taking derivatives with respect 
to t and evaluating at t = this becomes = F{u) + l'{u)R{u). Since l{u) is continuous on M^o 
with /(O) = 0, the above equation can be integrated to yield l{u) = — /J* ^||yds. By the first part of 
the proof this implies that a limit distribution L exists, with Laplace exponent l{u) coinciding with 
l{u). We conclude that also the probability laws L and L on Rj,o coincide, i.e. L = L . Conversely, 
assume that a limit distribution L exists. To show that L is also invariant, note that (|2.ip . ()2.6p 
and dZS]) imply 

Ptfu{x)dL{x) = exp u) - l{i^{t, u))) = exp I - / F(^(r, u))dr + / -^ds 



[0,00) y io ^0 

= { r = / fnix)dL{x). 

\Jo J J[0,oo) 

This completes the proof. □ 

Proof of Corollary \2.8[ Assume that Jj-j^ log^?Ti(d^) < 00. From the concavity of R{u), 
Lemma 12.41 and the fact that F{u) > for all u > we obtain that 

(A.4) 0<- r^ds<--^ r^ds = --^(bu+ r [ ^~^~'^ m{dOds] . 
^ ' ~ Jo R{s) - R'{0) Jo s R'{0) \ Jo 7(0,00) s ^ J 

In order to show that this upper bound is finite, it is enough to show that the double integral on 
the right takes a finite value. Since the integrand is positive, the integrals can be exchanged by the 
Tonelli-Fubini theorem. Defining the function M(^) = ^—^-^ds, we can write 

f f '^ m{dOds= [ M{Om{dO. 

Jo J{0,oo) ^ J(0,oo) 

An application of L'Hopital's formula reveals the following boundary behavior of M{(^): 

^ , M(e) , , M(£) 

A.5 lim — ^ = u and lim — = 1. 

^ log^ 

Choosing suitable constants Ci,C2 > we can bound M(^) from above by Ci^ on (0,1) and by 
C2 log.^ on [1, 00). Note that m(d^) integrates the function ^ i-^- Ci^ on (0, 1) by Theorem 12. !( and 
integrates the function 1-^ C2log.^ on [l,oo) by assumption. Hence 

[ M(C)m(de) <Ci [ e"i(dO + C2 [ log^m(dC) < 00, 

J{0,oo) J{OA) -^[l.oo) 

and we have shown that the upper bound in (1A.4P is finite and that (12. 7p holds true. 
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Suppose now that /^^-^ log ^ 77i(d.^) = oo. Since -R+(0) < we can find e,5 > such that 
R'_^iO) + e<0 and R{u) > iR'+{0) - e)u for all u G (0, 5). Hence, 

(A.6) - r^ds> r^ds = (bu+ r [ l^^m(de)ds i , 

for all u € (0, 5). Exchanging integrals by the Tonelli-Fubini theorem and using the function M(^) 
defined above we get 

f ^ m{d^)ds= [ M{i)m{di)>C2 f log ^ m(dO = cxd, 

(O.oo) 7(0,00) J[l,oo) 

where C2 > is a finite constant which exists by the second limit in ()A.5p . This shows that the 
right hand side of ()A.6P is infinite and hence that (j2.7p can not hold true. □ 

References 

Jean Bertoin. Levy processes. Cambridge University Press, 1996. 

Ma. Emilia Caballero, Jose Luis Perez Garmendia, and Geronimo Uribe Bravo. A Lamperti type 

representation of continuous-state branching processes with immigration. arXiv: 1012.2346, 2010. 
T. Chan, A. E. Kyprianou, and M. Savov. Smoothness of scale functions for spectrally negative 

Levy processes. Probability Theory and Related Fields, 150(3-4) :691-708, 2011. 
E. Qinlar and M. Pinsky. A stochastic integral in storage theory. Zeitschrift fiir Wahrscheinlichkeit- 

stheorie und verwandte Gebiete, 17:227-240, 1971. 
D. A. Dawson and Zenghu Li. Skew convolution semigroups and affine Markov processes. The 

Annals of Probability, 34(3):1103 - 1142, 2006. 
C. R. Heathcote. A branching process allowing immigration. Journal of the Royal Statistical Society 

B, 27(1):138~143, 1965. 

C. R. Heathcote. Corrections and comments on the paper 'A branching process allowing immigra- 
tion'. Journal of the Royal Statistical Society B, 28(1):213-217, 1966. 

Lancelot F. James, Bernard Roynette, and Mark Yor. Generalized gamma convolutions, Dirichlet 
means, Thorin measures, with explicit examples. Probability Surveys, 5:346-415, 2008. 

Z. J. Jurek and W. Vervaat. An integral representation for self-decomposable Banach space valued 
random variables. Zeitschrift fiir Wahrscheinlichkeitstheorie und verwandte Gebiete, 62:247-262, 
1983. 

Kiyoshi Kawazu and Shinzo Watanabe. Branching processes with immigration and related limit 
theorems. Theory of Probability and its Applications, XVI(l):36-54, 1971. 

Martin Keller-Ressel and Thomas Steiner. Yield curve shapes and the asymptotic short rate dis- 
tribution in affine one-factor models. Finance and Stochastics, 12(2) :149 - 172, 2008. 

A.E. Kyprianou and V. Rivero. Conjugate and complete scale functions for spectrally negative 
Levy processes. Electronic Journal of Probability, 13(57):1672-1701, 2008. 



LIMIT DISTRIBUTIONS OF CONTINUOUS-STATE BRANCHING PROCESSES WITH IMMIGRATION 19 

Andreas E. Kyprianou. Introductory Lectures on Fluctuations of Levy Processes with Applications. 
Springer, 2006. 

John Lamperti. Continuous-state branching processes. Bulletin of the AMS, 73:382-386, 1967. 
Zenghu Li. Measure- Valued Branching Markov Processes. Probabihty and its Applications. 
Springer, 2011. 

Ronnie Loeffen. On optimality of the barrier strategy in de Finetti's dividend problem for spectrally 

negative Levy processes. Annals of Applied Probability, 18(5):1669-1680, 2008. 
H. Masuda and N. Yoshida. Asymptotic expansion for Barndorff-Nielsen and Shephard's stochastic 

volatility model. Stochastic Processes and their Applications, 115:1167-1186, 2005. 
Mark Pinsky. Limit theorems for continuous state branching processes with immigration. Bulletin 

of the AMS, 78(2):242-244, 1972. 
Ken-Iti Sato. Levy processes and infinitely divisible distributions. Cambridge University Press, 

1999. 

Ken-iti Sato and M. Yamazato. Operator-selfdecomposable distributions as limit distributions of 
processes of Ornstein-Uhlenbeck type. Stochastic Processes and Applications, 17:73-100, 1984. 

M. Yamazato. Unimodality of infinitely divisible distribution functions of class L. Annals of 
Probability, 6:523-531, 1978. 

Department of Mathematics, TU Berlin, Germany 
E-mail address: mkeller@math.tu-berlin.de 

Department of Statistics, University of Warwick, UK 
E-mail address: a.mijatovic@warwick.ac.uk 



